Control of Multi-level Voltage States in a Hysteretic SQUID Ring-Resonator System 
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In this paper we study numerical solutions to the quasi-classical equations of motion for a SQUID 
ring-radio frequency (rf) resonator system in the regime where the ring is highly hysteretic. In line 
with experiment, we show that for a suitable choice of of ring circuit parameters the solutions to 
these equations of motion comprise sets of levels in the rf voltage-current dynamics of the coupled 
system. We further demonstrate that transitions, both up and down, between these levels can be 
controlled by voltage pulses applied to the system, thus opening up the possibility of high order 
(e.g. 10 state), multi-level logic and memory. 



I. INTRODUCTION 

In an earlier paper we reported on a new phenomena 
generated by the non-linear interaction of a SQUID ring 
(here, a single Josephson weak link enclosed by a thick su- 
perconducting ring) with a parallel LC resonant (tank) 
circuit. The block diagram for this coupled system is 
shown in figure ^ As has now become apparent, SQUID 
rings can display behaviour ranging from fully quantum 
mechanical through to quasi-classical [Sil&jiSiiSiI 

depending on the circuit parameters of the ring and the 
temperature (T) of the environment. At temperatures 
of a few K, and for relatively large weak link capaci- 
tances (10^12 - 10-"F), it is well established that 
a SQUID ring (inductance A, weak link critical current 
Ic) can be treated quasi-classically, i.e. as a particle mov- 
ing, with damping, in a cosine modulated parabolic po- 
tential. 
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FIG. 1: Block diagram for an inductively coupled SQUID 
ring-tank circuit system together with excitation, static flux 
bias and readout circuitry. 



where '^x is the external magnetic flux applied to the 
ring, $s is the total included flux in the ring and 
$0 ^— h/2e^ is the superconducting flux quantum. In 

this quasi-classical regime it is usual practice to de- 
scribe the SQUID ring, and its interactions, using the 
Resistivel y S hunted Junction plus Capacitance (RSJ-I-C) 
model U,m- Here, we denote the ring weak link effec- 
tive capacitance by Cs and the parallel link resistance by 
Rs. Furthermore, in making use of this model it is conve- 
nient to introduce the /3-parameter for the SQUID ring 

^= 27rA/c/$o^ ■ This allows us to distinguish between 
the parameter space in which Ic is always a single- valued 
function of (the inductive regime) and the hysteretic 
regime where Ic can be multi-valued in ^x 
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The cosine in the potential leads to a superconduct- 
ing screening current response in the SQUID ring which is 
manifestly a non-linear, and ^Q-periodic, function of ^s- 
When at least part of this external flux is time dependent, 
and applied through an inductively coupled resonant cir- 
cuit, the ring-resonator system displays non-linear dy- 
namical behaviour which has been the focus of much re- 
search over the last three decades 0, 0, 0| It is precisely 
this behaviour that forms the basis of the well know ac- 
biased SQUID magnetometer [1,E3E3- Although not 
invariably the case, it has been common practice to 
make use of a radio frequency (rf « 20MIIz) , parallel LC 
(tank circuit) resonator in these coupled magnetometer 
systems. In the work described here we adopt this fre- 
quency regime for the resonant circuit. Considered on its 
own, of course, the tank circuit circuit is strictly linear, 
i.e. on resonance the rf voltage {Vout) across it is linearly 
dependent on the level of rf drive input current (/„ (t)) 
applied. However, when a SQUID ring is coupled to such 
a tank circuit the situation can alter radically. With the 
SQUID ring potential above, and at finite critical cur- 
rent Ic, the screening supercurrent Ig = $s/A flowing in 
the ring to oppose any externally applied flux ^x is a 
non-linear function of this flux. Since both the SQUID 
ring and the tank circuit are macroscopic in nature, the 
back reaction through the inductive coupling leads to the 
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system displaying non-linear, even chaotic, dynamics. In 
general the stronger the non-linearity in i^x)i the more 
pronounced the non-linear behaviour can be. This non- 
linear behaviour is most easily seen in plots of Vout versus 
lin (t), these also being dependent on the level of static 
or quasi-static bias flux ($2,) applied to the SQUID ring. 
In SQUID magnetometer systems these characteristics 



display essentially constant rf voltage steps in the (time 
averaged) Vout versus (t) at intervals periodic in 
These SQUID steps are modulated in a <I>o-periodic man- 
ner by the external bias flux 's', 15, 1^. The ring-tank 
circuit system is described dynamically by two (coupled) 
equations of motion, one for the tank circuit and the 
other for the ring. These are given by 



+ i?;:^ + ^ + A(i-ir^) (^"""^ 



V + i?:ir + j + - (^QUID ring) (3) 



where the subscripts tc and s refer, respectively, to the 
tank circuit and the SQUID ring and fi (= M/Ltc) is the 
fraction of the flux coupled between the ring and the tank 
circuit. Thus, Ctc and Ltc are, respectively, the tank cir- 
cuit capacitance and inductance, is the flux in the 
tank circuit inductor, Rtc is the resistance of the parallel 

tank circuit on resonance and K {= ^ / KLt,^ quan- 
tifies the strength of the inductive coupling between the 
ring and tank circuit with a mutual inductance of M . 
In (|2Jl and Q the last terms on the right hand sides of 
these equations describe the back reaction between the 
ring and the tank circuit. 

In the original paper we used ^ and lO to com- 
pute the Vout versus dynamics of the ring-tank circuit 
system in the highly hysteretic, strongly underdamped, 
regime to model the observed experimental behaviour. 
What was recorded experimentally was a succession of 
plateau regions in the time averaged Vout versus lin char- 
acteristics of the ring-tank circuit system. Each plateau 
region consisted of a set a parallel steps at regular sep- 
arations along the lin axis. The lengths of these steps 
along this axis were much greater than those observed 
in standard hysteretic SQUID magnetometer character- 
istics, where, typically, /? ~ a few In addition, 
for each individual plateau the rf voltage {Vout) at which 
these steps occurred varied <I>o-periodically in <I>s. In ex- 
periment U the rin g-ta nk circuit system was observed to 
jump stochastically |l4ll5l[T6l| between the various steps 
associated with each plateau. Since, in practice, it was 
easy experimentally to access circuit parameters in which 
many steps per plateau region could be seen (a maxi- 
mum of 19 steps/plateau was recorded over the course of 
these experiments), there seemed a possibility that the 
plateaux and steps could be utilised to create multi-level 
logic as an alternative to the standard binary logic. We 
found that with simulated low temperature noise on the 
tank circuit drive current, introduced via a noise distri- 
bution for a thermal bath, the solutions to these coupled 
equations of motion modelled our experimental results 



very well, including the stochastic jumping between steps 
on particular plateaux. We argued in the paper P| that 
when the SQUID ring is sufficiently underdamped it no 
longer follows the potential adiabatically as the rf flux 
coupled in from the tank circuit changes with time. Es- 
sentially, it is this non-adiabatic response that generates 
the large multi-step plateaux in the ring-tank circuit Vout 
versus lin characteristics, each step corresponding to a 
different flux jump trajectory (local well to local well) in 
(^. Nevertheless, even though we had demonstrated that 
the underdamped RSJ-f^C description could model the 
experimentally observed plateaux and steps, the stochas- 
tic jumping between steps generated in the theoretical 
calculations required more explanation. In the experi- 
ments we took the view, we believed quite validly, that 
these jumping processes were caused by ambient noise. 
However, for technical reasons we were not able to mea- 
sure the full spectral density function of this noise and 
with the computer power then available to us we could 
not be certain that the stochastic jumping processes seen 
in the simulations did not arise because of computational 
inaccuracies. New calculations, at much higher accuracy, 
are the basis of the work reported here. At this new level 
we have been able to show that, from the viewpoint of 
the simulations presented in our previous paper T], the 
jumping processes arose due to the build up of compu- 
tational error in the numerical integration for the evo- 
lution of the system. Repeating the calculations, and 
eliminating this potential problem by using more accu- 
rate integration methods, only single steps are generated 
in each plateau region, even with large variance current 
noise added. Nevertheless, it was clear from experiment 
that the ring-tank circuit system could be perturbed to 
generate jumps. We now show how this can be achieved 
in a controlled way, pointing to possible device applica- 
tions in multi-level logic. 
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FIG. 2: Dynamical rf voltage (Vout) versus rf current (hn) 
characteristics for a highly hysteretic (/3 = 137) SQUID ring- 
tank circuit system with circuit parameters ring: Csq — 1 X 
10""F, Asq = 6 X 10"^°H, Rsq = ion, h = 75.2^A; tank cir- 
cuit: Ctc = 7.6x 10""'F, Ltc = 63nH and /i M^/LtcA;) = 
0.0087, a bare tank circuit resonant frequency 23MHz and 
= 0.0^0 (modulo n$o, n integer). 



II. SQUID RING-TANK CIRCUIT DYNAMICS 
IN THE HIGHLY HYSTERETIC REGIME 

With access to more computational power we consid- 
ered it important to establish the origin of the stochas- 
tic jumping, in the computer simulations as previously 
published Q. As an example of the possible origin of 
this jumping phenomenon we computed the Vout ver- 
sus lin characteristics for a SQUID ring-tank circuit 
in the large /3 regime. In this we chose suitable ring 
and tank circuit parameters for this regime, specifically 
for the SQUID ring: A, = 6 x lO'^^R, 1^ = 75.2/zA 
iP = 137), Cs = IQ-^^F, R, = ion and for the tank 
circuit: Ctc = 7.6 x 10~^"F, Ltc = 63nII and a quality 
factor Q of 500, the latter parameters yielding a bare 
(uncoupled to the ring) tank circuit resonant frequency 
of 23MHz. We also set (= M^/UcA; ) = 0.0087, a cou- 
pling quite typical of low noise SQUID ring-tank circuit 
systems. These are the circuit parameters used through- 
out this paper. In figure|21we show two sets of computed 
Vout versus characteristics using these ring and tank 



circuit parameters. With these values of Rs and Cs the 
SQUID is underdamped. The originally reported char- 
acteristic 0, calculated at relatively low accuracy using 
a fourth order Runge-Kutta numerical integration rou- 
tine with an adaptive step size algorithm, is plotted in 
figure |2Ia) for the first plateau region. In generating this 
characteristic 4.2 Kelvin current noise, with a thermal 
noise distribution, has been introduced and the bias flux 
has been set at $a: = 0.0$o (modulo n$o). Here, the 
multi-step solutions on the first plateau, and the stochas- 
tic jumping between these, are perfectly clear. Since the 
initial calculations were performed, we have had access 
to much greater computational power and this has al- 
lowed us to apply significantly higher accuracy to our 
numerical integration. Subsequently we have found that 
this greatly improved accuracy leads to the suppression 
of the jumps between the multiple levels in the Vout ver- 
sus lin characteristics. This is very apparent in the other 
solution shown in figure [2Ib) for which we used exactly 
the same SQUID and tank circuit parameters as in the 
characteristic of figure[2Ia). This second solution implies 
that the experimentally observed jumps may be driven 
by noise processes but not through current noise as we 
previously thought. Within the circuit model we have 
adopted, the alternative is that the driving force arises 
from some form of voltage perturbation. In order to in- 
vestigate this possibility we present a number of simu- 
lations of the response of this system to voltage pulses 
applied to the SQUID ring. 

In considering the dynamical basis for the plateaux 
(and the step levels) it is helpful to examine the way in 
which the SQUID ring flux (<I>s) changes in time with 
the rf drive current (/i„). Such a response, is shown 
in figure 01 for a maximum (sinusoidal) drive flux in the 
tank circuit of 13.5<I'o peak to peak (equivalent to a drive 
flux in the SQUID ring of 0.1175$o) ■ In this figure the 
SQUID ring flux is normalised to $o and the time axis 
is plotted in units of the reciprocal tank circuit period 
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FIG. 3: Time dependence of the normalised SQUID ring flux 
^sq/^o as a function of time, normalised to the tank circuit 
period LtcCtc- 
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FIG. 4: The SQUID ring potential for the system is shown in grey for (a) a positive $3 going, outward, muIti-4>o trajectory; 
(b) return trajectory for (a); (c) a negative going, multi-$o trajectory; (d) return trajectory for (c). 



(2TTy/LtcCtc) and, for comparative timing purposes, we 
show the form of Im in Hght grey. What this computed 
plot demonstrates is that, for the most part, the oscil- 
lating flux (<l>s) in the underdamped SQUID ring corre- 
sponds to the ring being confined to a single well. In fig- 
ure |3| this single well SQUID ring flux response is shown 
oscillating over small amphtudes around $s/^'o ~ 0. 
However, when the drive current reaches maximum am- 
plitude it is found that the SQUID ring can move outside 
the confines of a single well and make more extended tra- 
jectories in the potential i.e. the flux in the ring can 
jump non-locally (and non-adiabatically) several $0 in 
^s-space between initial and final (target) wells 0|. We 
now provide some illustrative examples of stable solutions 
of equations 10) and Q which give rise to the different 
voltage levels seen in figure l^^a). 

In the particular example of figure 13 these non-local 
jumps take place over ± six wells in $s-space relative 
to the lowest energy well in the potential. Transposed 
to the rf dynamics of the coupled system, each multi-$o 
traversal leads to one of the set of plateau step levels in 
Vout versus lin , as exemplified in figure |21 Examples of 



these multi-$o trajectories (jumps) in the SQUID ring 
potential can be seen in figure ^ again for the system 
of figure m Taking the origin as the lowest well in this 
potential, the SQUID flux trajectories shown in figures|2| 
and ^correspond to the level denoted Iq in figure|2Ia). If 
the origin is displaced to the next lowest well in the po- 
tential, the Vout versus lin characteristic of the coupled 
system moves down by one level from the original solu- 
tion. This level is denoted li in the solutions plotted in 
figure Ufa). Moving the origin again to the next lowest 
well repeats the process, i.e. the dynamics now corre- 
spond to the level denoted by I2 in the solutions plotted 
in flgure|2Ia). Again, moving the origin to the next lowest 
well generates the level denoted I3 in figure Efa). How- 
ever, in our example of a 6$o excursion this level {I3), 
is equivalent to having displaced our origin to the mid- 
dle well. We find that beyond this central well in the 
excursion, attempting to localise our trajectory around 
one of the remaining three wells simply reverses the shift 
between the levels on the plateau in flgure|2Ia), i.e. start- 
ing at I3 and returning to Iq through levels I2 and li. To 
illuminate this discussion, we can examine, for level Iq, 
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the excursions in the potential shown in the figure 0] 
Starting in the absolute minimum well of the potential 
at $s/$o = 0, with = (modulo n$o), the ring can 
either execute local motion in this well or occasionally 
move several $0 to a target well. With the ring potential 
plotted in light grey in figureQ] we show by computation 
in figure Efa) one such extended trajectory, starting in 
the lowest well and followed by a positive going trajec- 
tory in $s. In this example, having reached the target 
well the ring completes almost half a tank circuit period 
in this well before returning to its original starting point 
in $s, as shown in figure ^b). This mirrors the dynam- 
ical behaviour shown in the first four jumps of figure |3| 
It is also important to note that these non-local excur- 
sions do not occur every time the rf drive flux reaches its 
maximum amplitude, although as the amplitude of the rf 
drive grows this process happens more frequently. Fur- 
ther multi-$o excursions in the ring potential, but now 
for negative going can be seen in figures QJ^c) and (d) 
for the outward and return paths, respectively. These 
correspond to the last two pulses in figure |2| In general 
the target well in the SQUID ring is not symmetric about 
its centre (minimum) and this is seen to affect the life- 
time of the ring in this target well, as evidenced in plots 
presented in figure 01 In essence the consequence of this 
is that the ring can reach a new traversal (target) point 
in drive flux before it has completed half a $0 period in 
^s- This phenomenon, arising from the localisation in a 
non-symmetric well in the SQUID ring potential, means 
that the excursions in <I>s are only made in one direction. 

With the details of figure 01 in mind, and as a possi- 
ble means to navigate the solutions of 10) and ) , we 
now consider the effect of rapidly changing voltages ap- 
plied to the SQUID ring. In order make this relatively 
simple, and physically transparent, we simulate the ap- 
plication of appropriately shaped voltage pulses to the 
SQUID ring, and follow the tank circuit voltage (Vout) 
response. In our example we shall now consider positive 
amplitude trapezoidal voltage pulses (inset in figure \^ 
with an upper voltage state time duration of 0.01 tank 
circuit periods which is activated during the flux traver- 
sal shown in bold in figure 01d). We found that it was 
only within such highlighted activation regions that the 
system could be made to respond to these voltage pulses, 
leading to a level change in Vout- For reference we term 
pulses with positive amplitude as A type pulse and, con- 
versely, negative amplitude pulses are denoted as B type. 
In figure [3 we have plotted the computed tank circuit 
voltage response as a function of normalised tank circuit 
time (i.e. in terms of the tank circuit period) for the 
rf drive amplitude (/i„=222nA) set in the middle of the 
first plateau for the level denoted Iq in figureEIa). In this 
calculation there was no extraneous flux or voltage noise 
present and the system was allowed to reach a steady 
state within this plateau (we note that the initial sharp 
rise at the beginning of the voltage dynamics is due to 
transient behaviour). However, after a set time, denoted 
Ai in figure [S] a voltage pulse of the form shown inset 
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FIG. 5: Vout versus normalised time characteristics for the 
SQUID ring-tank circuit system of figure |21 showing the con- 
trolled step to step jumping which can be induced using suit- 
ably configured positive (A-type) ramp pulses in the absence 
of extraneous voltage noise; here we show the effect of in- 
creasing the ramp amplitude to induce jumps over one, two 
or three step intervals. 



in this figure was applied. Here, the amplitudes of the 
pulses Ai,i = 1, 2, 3 used were 15, 75 and 120^V, respec- 
tively. As is demonstrated in figure [S] the application of 
a single voltage pulse can change the voltage level of our 
system by one or more levels at a time, as summarised 
in table m If the noiseless (or, in practice, a low enough 



Pulse Type 


Lower Bound 


Upper Bound 


Value Used 


Ai 


15 


55 


15 


A2 


56 


90 


75 


A3 


91 


126 


120 



TABLE I: Summary of pulse amplitude ranges (/^V) as de- 
scribed in the text. 

voltage noise) situation can be realised experimentally, 
the computational results of figure |31 indicate that multi- 
level logic, based on SQUID ring-resonator systems could 
be a feasible proposition. If so, it is reasonable to assume 
that other forms of voltage pulse could be used for this 
purposes. With regard to our example, we note that the 
number of levels on the plateaux depends on the ring 
parameter values. In the example of figure |S1 this is four 
but could be significantly larger. To further illustrate the 
control possible we show in figure El the result of apply- 
ing two A-type pulses, in time sequence, to the SQUID 
ring-tank circuit system. 

Of course, the activation pulse in voltage can be nega- 
tive (B-type) as well as positive. The effect of using such 
a pulse after an A-type pulse is shown in figure with 
the form of the B-type pulse shown explicitly in the inset 
of this figure. This shows that in the absence of extrane- 
ous noise, and by a suitable choice of voltage pulse form, 
we can move at will between the steps on any particular 
plateau in Vout versus To demonstrate that we can 
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FIG. 6: Vout versus normalised time ciiaracteristics for the 
SQUID ring-tanlt circuit system of figure El showing the effect 
of two sequential A-type voltage pulses in the absence of ex- 
traneous voltage noise; as is apparent these step jumps can 
be induced in a controlled manner. 
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FIG. 8: Vout versus normalised time characteristics for the 
SQUID ring-tank circuit system of figure |21 in the absence of 
extraneous voltage noise showing the efi'ect of applying se- 
quentially at difi^ering time intervals as a set A-type voltage 
ramps. 
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FIG. 7: Vout versus normalised time characteristics for the 
SQUID ring-tank circuit system of figure|5|in the absence of 
extraneous voltage noise showing the effect of applying se- 
quential A and B-type voltage ramps. 




FIG. 9: Graphical representation in the manner in which the 
wells in the SQUID potential correspond to voltage levels in 
the Vout versus Ii„ characteristic of a highly hysteretic SQUID 
ring-tank circuit system. The vertices around the polygons 
indicate which well (in a six well example - see figure 0]l the 
SQUID oscillation is localised in, whilst the letters on the 
projection denote which one of the set of levels on a particular 
plateau in Vout versus /i„ the coupled system occupies. See 
text for more details. 



induce step to step jumping at essentially any time of 
our choosing, we show in figure|Sla sequence of jumps in- 
duced by a sequence of A-type pulses spaced by different 
time intervals. We note that the application of the last 
pulse causes the tank circuit voltage to step up rather 
than down. 

The computed solutions of figures [3 to |H1 show very 
well the level to level jumping induced by various pulse 
sequences applied to a highly hysteretic (large (3) SQUID 
ring-tank circuit system. These solutions, which are the 
end result of very non-linear interactions between the ring 
and the tank circuit, may be appreciated more clearly by 
the following qualitative argument. As regards the jump- 
ing process between levels on any particular plateau in 
Vout versus Ii„, the system operates cyclically from an 
initial (local) well in the SQUID potential (for example 
the wq well in the potential of figure^. This is shown di- 



agrammatically in figure |5| The way in which movement 
between wells in the SQUID potential relates to a given 
voltage level within a particular plateau can be deter- 
mined by thinking of each well being mapped on to one 
of the vertices of a polygon (where the order of the poly- 
gon is dependent on the parameter values of the ring-tank 
circuit system under investigation). This is illustrated in 
figureinia) for the system considered throughout this pa- 
per with the vertices labelled Vi for i = 0, 5 and where 
the mapping is applied cyclically from the lowest (ref- 
erence) well of the potential on to each vertex (for our 
example using the labels for the wells Wi in figureEJa) we 
have the mapping Wi "fimode for i G Z). Here, we can 
find the particular tank circuit voltage level associated 
with each vertex by projection on to the horizontal axis. 
For example, the vertices uiand U5 correspond to having 
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the SQUID ring-tank circuit system in level h in figure|21 
Clearly, therefore, if the system is moved through an ap- 
propriate number of wells it will return to the same state 
of operation in which it began. 

Using this graphical map as the guide, in this paper we 
have investigated a system with 4 levels, corresponding 
to figure IHl^a). We note that in this example the system 
must traverse 6 wells in the potential before it returns 
to its original state of operation. If we choose differ- 
ent SQUID parameters the system can access a variable 
number of wells within its cycle of operation. However, 
the system will always display a symmetry about the 
horizontal axis. For comparison, we show in figure |5Jb) 
a diagram illustrating the situation where the ring-tank 
circuit system traverses an odd number of wells in the 
SQUID potential, in this case five. 

By the application of suitable voltage pulses the sys- 
tem can be moved around the vertices of the polygon, 
so translating the SQUID ring from one potential well to 
another. Thus, in figure EJa) apphcation of positive (A- 
type) pulses causes the ring to move around the polygon 
in a clockwise direction, whilst negative (B-type) pulses 
lead to anti-clockwise movement. This movement around 
the polygon can be performed utilising nearest neigh- 
bour, single well, translations, or by multi-<I>o traversals 
between non-local wells in the potential, depending on 
the magnitude of the applied pulse The actual level in 
Vout versus /i„ into which the system settles is deter- 
mined by three factors: (i) the well in which the SQUID 
ring is currently localised (ii) the size of voltage pulse 
used to change location in the potential and (iii) the di- 
rection of traversal of the polygon (whether A or B-type 
pulses are utihsed). 

III. CONCLUSIONS 

In this work we have demonstrated that, in principle, it 
is possible to use highly hysteretic (large (3) SQUID ring- 



tank circuit systems as the basis for multi-level logic or 
memory devices Following earlier ex- 

periments, where giant SQUID magnetometer plateaux 
were observed containing sets of constant voltage (Vout) 
steps, we have shown that the jumping processes between 
steps on a particular plateau can be generated by voltage 
pulses applied to the system. In the absence of noise (in 
our computational modelling achieved by the use of high 
accuracy numerical integration techniques) , we find that 
a high /3 SQUID ring-tank circuit system can remain sta- 
bly on one of the steps on any particular plateau until 
a suitably shaped voltage pulse is applied to the ring. 
Following such a pulse, the ring can be made to jump 
in a controlled manner to other steps at different lev- 
els of Vout- This new level can be either above or be- 
low the original voltage level, depending on the present 
state of the system and the type of pulse used. Pro- 
vided it is possible to reduce the ambient voltage noise 
on the SQUID ring sufficiently, for example by the use of 
cryogenically cooled GaAsFET or HEMT-based pre- 
amplifier electronics |23|, it should prove possible to de- 
velop devices with a controlled voltage response which 
can be selected, and modified, at will. This could prove 
useful in such areas as multi-level logic or finite state ma- 
chinery [Illlllllllllll. As part of any such develop- 
ments we would expect to see interesting new fields open 
up in the non-linear dynamics of these highly hysteretic 
systems. 
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